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THE PROCESS OF SEPARATION IN THE TURBULENT. FRICTION LAYER* 

By E. Gruschwitz 



Tlie separation of the flow which occurs at large an- 
gles of attack on the suction side of an airplane wing is 
attributable to phenomena in the flowing fluid layer ad- 
jacent to the surface; the fluid particles slowed up by 
the friction on the surface can no longer advance against 
an unduly great pressure rise. It is of vital importance 
that there exist two types of flow - laminar and turbu- 
lent - in the fluid layer flowing in the immediate vicin- 
ity of a body. According to Prandtl, by whom the whole 
theory was developed, we speak in the first case of a lam- 
inar boundary layer, in the second, of a turbulent fric- 
tion layer. 

Beginning at the stagnation point, there is always a 
laminar boundary layer present, which at some distance 
from the stagnation point changes into a turbulent fric- 
tion layer on sufficiently large bodies and by sufficiently 
high speeds. In view of the overcoming of a rise in pres- 
sure, the turbulent flow attitude is markedly more favor- 
able than the laminar attitude, since the turbulent mix- 
ing effects a much more intensified impulse transfer on 
the decelerated fluid particles than the mere viscosity. 



I. METHOD 3? OR APPROXIMATE CALCULATION OP FRICTION LAYER 



In another report** the writer attacked the problem 
of making the turbulent friction layer amenable to cal- 
culation for the case of steady plane- flow. Being prac- 
tically ignorant of the nature of the turbulent isolated 
processes it, of course, pertained only to - an approxima- 
tion method derived from experiments-. ; The method is brief- 
ly, as follows: 



" it 

*"Uber den Ablo sungsvorgang in der turbulenten Reioungs- 
schicht." Z.F.M., June 14, 1932, pp. 308-312. 
**Ingenieur-Archiv II (1931), p. 321. 
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The premises of all 'boundary-layer calculations are 
that the speed at relatively short distance from the sur- 
face becomes" the speed of a potential flow and that within 
the friction layer the static pressure perpendicular to 
the -surface does not change. Admittedly, only the mean 
values with respect to time in the turbulent flow are of 
interest so long as it is not desired to follow the tur- 
bulent isolated motions. The velocity distribution in 
the friction layer past a flat plate exposed longitudinal- 
ly to a flow, that is, in the absence of eny pressure gra- 
dient, can be interpolated with close approximation by 
power formula. (Reference 1.) lor velocity distributions 
which reveal the friction layers in the presence of great- 
er pressure gradients along the body, such interpolation 
is impossible, as Figure 1 shows. Besides, it appears 
hopeless of finding a practical mathematical term for a 
set of curves whose individual curves would reproduce the 
velocity distributions. One therefore must forego such 
an interpolation.* 

One first legitimacy., which must suffice the friction 
layer, is supplied by the momentum theory, as first applied 
by v. Earman to boundary layers. The derived equation reads: 



dx 2 q dx P tf 3 K 1 

where 

00 (W - w)w l r °° , w\ , 

* = I — dy and H = * I K " *J dy 

(x = length of arc at surface of body in direction of flow, 
y = coordinate at right angle to it, W = magnitude of po- 

tential flow, q = i^-— » w = speed in friction layer, t 0 = 

shear stress at wall and P = density). 

*In a recently published article (Horst Muller, "IFricti onal 
Resistance of Bodies ill a Flow," '7erft, Reederei, Hafen, 
1932, Ho, 4-, p. 54), the velocity profiles in friction lay- 
ers by pressure rise and drop are approximated by the for- 
mula i = 0p; /? (6 = distance in which the friction layer 

changes to potential flew, and coacln^i ons drawn (after 
the pattern 'of v, Karman, loc* cit.) for the friction lay- 
er on a flat plate. It yields a practical estimate, of the 
frictional resistance in comparatively simple manner. 
Since a solid profile form is assumed, it yields, of course, 
no way to ihe separation problem. 
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The quantity £, having', the dimension of a length, 
can he* looked upon as .a criterion for the thickness of the 
friction layer; # P W 3 is the loss of momentum due to 
friction in the friction layer. In the absence of this 
layer, that is, in the presence of the potential velocity 
W directly on the surface, the flow of momentum between 
the surface and a point at distance $ in unit of time 
would "be as great as the loss through friction in the 
friction layer. 

Equation (l) is essentially a differential equation 
for £, with, however, H • and T 0 /PW 3 as further un- 
known quantities. Compared to the other terms of the 
equation, on the other hand, quantity T 0 /PW 2 is quite 
small when the pressure gradient 




assumes an appreciable value. Consequently, it suffices 
to substitute a constant approximation value. It was 
found that T 0 /pW 2 = 0.002 was quite acceptable.' Quanti- 
ty E shall be referred to again later. 

It being important to investigate the processes which 
lead to separation of the friction layer, it does not suf- 
fice to follow up its thickness, but rather to observe the 
form of the velocity profiles. We concentrate on the ve- 
locity at a point of the friction layer lying near the 
wall. When this velocity drops below a certain value as 
compared to the velocity outside of the friction layer, 
it may be conceded that separation is imminent. To define 
this point we are without all means except that selected 
for the thickness of the friction layer and, in particu- 
lar, it has proved advisable to choose the point at dis- 
tance # from the wall. & is quite small compared to the 
width taken up by the friction layer; for example, by a 

velocity distribution of £ = (''£V /7 £ ' - ;~r 6 (the ap- 

proximation formula for the flat plate). 

The fluid particle at this point, whose velocity may 
be designated at w lf can be expressed by an approximate 
motion equation. The turbulent motions transfer impulses 
from the outside to the particles and give off impulses 
on the fluid layer farther inside. Of these impulses it 
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is assumed that they are determined "by w lf ¥, £ and den- 
sity P (possibly the viscosity might he of influence, but 
in turbulent motions it is mostly negligible compared to 
the turbulent momentum exchange). Furthermore, it should 
be borne in mind that the curve y = on which the 

fluid particle is visualized as moving, has practically 
the aspect of a streamline. Then the equation of motion 
reads : 

p w * + dl = f < w i» w ' *• p > ( 2 ) 

PW3. 2 

(p = statio pressure) or, putting — - — * p = g x 

|§i = f (w lt W, 4, P) (3) 



Since every physical relation must lend itself to being 
expressed as equation between diinen sionle s s quantities, 
and inasm uch as i&r the five quantities appearing in (3) 
(on ly t wQ^ permit the formation of ^independent dimensionless 
qua!TET : ries7~ equation (3) can also be written in the form 




■ * (r) M 



Measurements of the velocity distribution in fric- 
tion layers by pressure rise and pressure drop have actu- 
ally revealed in close approximation a relationship of 
form (4), which may be interpolated by the equation 



* dg 



l 



q d x 



= 0,00394 T\ - 0.00461 (5) 



whereby - /Wj.\ a 



_ _ n /"IV 



(6) 



(Fig. 2.) It is seen from (5) that the turbulent motions 
extract more momentum from the fluid particle with veloc- 
ity W]_ than it obtains from without when 

Tl < 0^00461 = 0 51g 
1 0.00994 U « UJ - D ' 

whereas the momentum from without preponderates when T) > 
0.516, that is, so much more as r\ is greater, or in 
other words, as the particles remained behind relative, to 
the potential flow. Equation (6) can also be written: 
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or 




Since q + p is a constant quantity, it follows: 

agi = a(q'n). 
d x d x ' 

I 

i 

so that (5) can "be written as 

A iLLaZll = - 0.00394 q T| + 0.00461 q (7) 
d x 

Visualize w/T7 plotted against y/# and assume 
that the thus produced sheaf of curves is approximately 
of one parameter, so that quantity w x /W or T) in its 
stead, can fee considered as sheaf parameter. Then the 
quantity 

h - \f ( i - ^* , - / (i - f) a (I) 

occurring in (1-) must be a function of T) , as is actually 
the case according to the experiments. (See fig. 3.) Be- 
tween H, S- and the length 

sometimes used as criterion for the thickness of the fric- 
tion layer, there exists moreover, the relationship H = 
5 a /$, as seen from the definition of H. 

The result is a system of differential equations in 
(l) and (7) from which "by known velocity past a body and 
by known initial values the friction layer can be comput- 
ed, because T) defines the form and & and the prescribed 
potential velocity W, the expanse of the velocity pro- 
file. The presumption hereby is that the thickness of the 
friction lay-ex .is small pompared to the curvature radius 
of the surface, because (7) has been derived from measure- 
ments on a flat wall, whereby the rise and drop in pres- 
sure had been produced by corresponding formation of an 
opposite wall.. On markedly curved walls a certain curva- 
ture effect .'is to be expected according to the experiments 
which are still being continued. As concerns,- the initial 
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values, they stipulate an idealization, if the turbulent 
friction layer is to "be computed from the beginning. It 
is assumed that the change from laminar "boundary layer 
into turbulent friction layer occurs with a sudden jump. 
In reality, of course, this change-over extends over a 
great distance, but it has been found that this idealiza- 
tion is very expedient when placing the point of transi- 
tion where the boundary layer begins to become turbulent, 
and there to give the turbulent friction layer the value 
r\ = 0.1. As initial value for one must select the 

value which the laminar boundary layer exhibits at its 
end, to ensure steady impulse content. 

A close approximation for # can be obtained by sub- 
stituting constants (say, 

. H = 1.5, -X- = 0.002). 

P W 

T 

for H and — in equation (l) which, thus becoming a 
p TT 

linear differential equation of the first order for this 
approximation, is inserted in (7) to yield a linear dif- 
ferential equation of the first order for q T). After 
computing this, the course of H can be ascertained and 
written in (l) , which yields a linear differential equa- 
tion of the first order for a second approximation of # 
for computing a second approximation of T), etc. The res- 
olution of the linear differential equations is best done 
by E. Czuber's graphical method. (Reference 2.) In the 
examples computed heretofore, it was never necessary to go 
beyond the second approximation. In many cases even the 
first showed satisfactory agreement with the experiments. 
Separation is expected when T) > 0,80. 

The objectional feature of this method is that it 
stipulates the knowledge of the laminar boundary layer 
and of the point of transition into turbulence. For ex- 
ample, having " comput ed the flow around an airplane wing 
according to the potential theory, it is necessary to com- 
pute the laminar boundary layer, which can be effected 
approximately by the Karman-Pohlhausen method. (Refer- 
ence 3.) The point of transition is primarily contingent 
upon the value of R = WS-fv, that is, on a Reynolds Slum- 
ber formed with the quantities decisive for the expansion 
of the velocity profile. However, no exact figure of S, 
at which the boundary layer becomes turbulent, can be 
given. Obviously, it also depends on the form of the lam- 
inar velocity profiles and undoubtedly also on the prop- 
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erty of the air stream to which the airplane wing is ex- 
posed,' According—to measurements hitherto, the laminar 
boundary layer obtains up to E = S50. If return flow 
occurs before reaching this figure, the transition to tiir- 
bulence occurs at the point of inception of this return 
flow, provided R has there exceeded a certain figure 
which, although not definitely determined, should be below 
250, since that is the lowest figure for R observed in 
all turbulent friction layers up to now. If the boundary 
layer evinces no, sufficiently high. E at the point of in- 
cipient return flow, there is no formation of turbulent 
friction layer and separation occurs. 

II. INFLUENCE OF REYNOLDS NUMBER Oil' EE .!. C TION LATER 

It is remarkable that the Reynolds Number R = TV &/v 
does not occur in equation (7) . For its derivation we em- 
ployed measurements in friction layers in which R = 250 
to 5,450 occurred, Whether or not there is a relation- 
ship of R in experiments on friction layers with higher 
Reynolds Numbers, is an open question, though improbable 
because of the already mentioned subordinate role which 
viscosity plays in comparison with the turbulent motions. 
In (l) the influence of the Reynolds Number is confined 
to the value of t 0 /PF 2 . But this influence is quite 
small, and T 0 /p¥ 2 is very small compared with the other 
two terms when dp/dx is great, as i s the case in flows 
in which separation i s in question, so that as a result, 
the dependence of this equa-tion on R can be ignored in 
such flows, and the turbulent friction layer can be con- 
sidered as being independent of- R. 

Nevertheless, the total process is not independent of 

W t 

the Reynolds Number R% = — which, as customary;, is 

formed with the flow velocity W 0 . . and a characteristic 
length t of the exposed body (as' with the chord, say, on 
the airplane wing). The laminar boundary layer, always 
present at the beginning is, of course, dependent On. the 
via c»o. si ty . a n A.j consequently, on any formulated Reynolds 

Number; specifically, i5-/t is proportional to j=~? (ref- 

v Rt 

erence 4) , according to boundary layer theory. The rela- 
tive thickness of the boundary layer at a certain point 

decreases as Rt increases, whereas the value of R = 
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increases, at this, p bint when Et increases., tie cause &/t 

proportional to ' -~rr means that . E ; .= -is proportion- 

. ;. • ,'. v • Et , ' it _ 

al to — ." ■ ' - that is, proportional to V E t* According- 

ly, as' Et increases, the point of transition to turbu- 
lence moves forward, while on the contrary, the friction 
layer begins with a smaller # value.-. .Equation (5). re- 
veals that the. impulse given a particle with, velocity w x , 
positive when r\ > 0,516, and thus aiding to overcome 
the pressure rise, is inversely proportional to so 
that a smaller & value. is more favorable in the sense of 
avoidance of separation (r\ > 0,516, ordinarily by pres- 
sure rise). The shifting of the. transition point, act s in 
the reverse sense. The farther forward the beginning of 
the turbulent friction layer, the greater the pressure 
rise it has to. overcome and the . more ^mminent the danger 
of separation, .,' 

In this way one can visualize . the manner in which the 
dependence of the maximum lift of an airplane wing on the 
Eeynolds Number- Et comes into being. It may be assumed 
that the maximum lift occurs when; the separation begins at 
the trailing- edge or is- directly imminent there.. In wings 
"with well-rounded- of f leading edges and fairly thick' pro- 
file, the turbulent transition occurs considerably behind 
the pressure minimum at lower Eeynolds Numbers. For high- 
er Eeynolds Slumbers the transition point shifts forward" 
and the flow separates therefore at lower-angles of attack, 
i.e., the coefficient of "maximum lift becomes smaller. 
Lastly, by further increase in Eeynolds Number the point 
of transition finally falls in the vicinity of the pres- 
sure minimum, wher# it probably remains or at least shifts 
very little, because in the zone ahead of it, the bound- 
ary layer will be so thin by pronounced pressure drop as 
to scarcely induce turbulence even by very high Eeynolds 
Numbers. Then the influence of the diminution on £ would 
have to make itself felt, and an increase in maximum lift 
is surmised. In wings" with sharp leading edge, that is, 
such with thin profile, especially, the leading edge evi-.. 
dently produces the transition to turbulence, and since-- ;• 
it then always remains at the same place, .only the dimi- 
nution of $ is of influence by an increase' in Eeynolds 
Number, so that a slight increase in maximum ..lift ^coeffi:- ■ 
cient is. to be expected. As a matted of. fact, comparative 
measurements at different Eeynolds Numbers such as those 
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made at Gottingen* and in the variable density wind tunne*H 
of ..the National Advisory Commit t.e.e for Aeronautics (ref- 
erence 5), reveal that the maximum lift coefficient ordi- 
narily decreases "by increasing Reynolds Number for thick 
profiles with well.-r ounded-of : f leading edges',, and a slight 
increase for thin profiles. But for uno.b j ec ti onable proof 
of the surmise, according tp t which the coefficient of max- 
imum lift should increase -.at high Reynold's Numbers even 
for thick profiles, the 'measurement s made; heretofore are 
insufficient. 

A very similar process occurs on cylinder a. and spheres 
also. Admittedly, equation's (l) and (7) are invalid for 
the sphere, in so far as it pertains to a symmetrical ro- 
tational flow. But qualitatively the expected behavior of 
the friction layer is similar. Resistance measurements .on 
cylinders and spheres reveal that beginning from the crit- 
ical value their resistance coefficient increases. (Ref- 
erence 6.) By "critical value" is meant the value of the 
Reynolds Number W 0 d/u (¥ 0 = flow velocity, d = diame- 
ter) at which a tu.rbulent friction layer occurs. The rea- 
son for this is, according to the above, the forward dis- 
placement of the !• trail si tion -point and with it its point 
of separation by increasing Reynolds Number, with a re- 
sultant enlargement of eddy zone and consequently, resist- 
ance • 



Translation by J, Vanier, 
National Advisory Committee 
for Aeronautics., 



*See Brgebnisse der Aer odynami schen Ver suchsanstalt zu 
Gottingen, Report I, p f 54. 
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y, mm 

Fig. 1 Velocity profiles on G-6'ttingen airfoil 
No. 387 (chord = 40 cm), flow velocity 
W 0 = 30.8 m/s, angle of attack a = 12°, x = 
length of arc along suction side from foremost 
point (total length to trailing edge = 42.7 cm), 
y = distance from surface, w = velocity in 'bound- 
ary layer or friction layer. The first two ve- 
locity profiles are laminar, the rest, turbulent. 
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